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CLT for the zeros of Kostlan Shub Smale random polynomials 


Federico Dalmao*^ 
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Abstract 

In this paper we find the asymptotic main term of the variance of the number of roots of Kostlan- 
Shub-Smale random polynomials and prove a central limit theorem for the number of roots as the 
degree goes to infinity. 


Resume 

Dans ce papier nous trouvons le terme asymptotique dominant de la variance du nombre de racines 
reels des polynomes aleatoires de Kostlan-Shub-Smale et demontrons un theoreme de la limite 
centrale pour ce nombre de racines. 

Consider the Kostlan-Shub-Smale (KSS for short) ensemble of random polynomials: 

d 

Xd{x) := X ^ X e M, 

n=0 

where d is the degree of the polynomial and the coefficients (alfb independent centered Gaussian 
random variables whose variances are the binomial coefficients, more precisely Var(alfb = (n)- 
Denote by the number of real roots of X^, that is 

Nd := #{x e M : Xd{x) = 0}. 

It is well known that K{Nd) = y/d [10, 15]. The aim of this paper is to prove the following result. 

Theorem 1. The variance of the number of real roots Nd of KSS random polynomials verifies 

Var(Arf) 2 

hm - -j= — = cr , 

d-^OQ y/d 

with 0 < < oo given in Proposition 1 ~ 0.57 ...). Furthermore, Nd verifies the CLT 


Nd-Vd 

di/4 


^ A(0;cjh. 

d—)-oo 


The number of roots of random polynomials has been under the attention of physicists and math¬ 
ematicians for a long time. The first results for particular choices of the coefficients are due to Bloch 
and Polya [5] in 1932. After many successive improvements and generalisations, in 1974 Maslova [12] 
stated the CLT for the number of zeros for polynomials with i.i.d. centered coefficients with finite vari¬ 
ance. For related results see the review by Bharucha-Reid and Sambandham [6] or the Introduction 
in [2] and references therein. 
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The study of Kostlan-Shub-Smale (m x m systems) of polynomials started in the early nineties 
by Kostlan [10], Bogomolny, Bohigas and Loboeuf [4] and Shub and Smale [15]. The mean number of 
roots [15, 10], some asymptotics as m —)■ oo for the variance [16] and for the probability of not having 
any zeros on intervals [14] are known. See also the review by Kostlan [11] and references therein. 

We restrict our attention to the case m = 1. The mean number of real roots is y/d. This fact 
shows a remarkable difference with the polynomials with i.i.d. centered coefficients which asymptotic 
mean number of roots is 21og((i)/7r. 

Our tools are the Rice formulas for the (factorial) moments of the number of roots [3, 7]; Kratz- 
Leon’s version of the chaotic expansions for the number of zeros [8], Kratz-Leon method [9] and the 
Fourth Moment Theorem [13]. This method has been applied to trigonometric polynomials by Azais, 
Dalmao and Leon [2]. 

The paper is organised as follows. Section 1 contains some preliminaries and sets the problem in a 
more convenient way. Section 2 deals with the asymptotic behaviour of the variance of N^. In Section 
3 the asymptotic normality of the standardised is obtained. Section 4 contains the proofs of some 
auxiliary lemmas. 

1 Preliminaries 

By the binomial theorem, the covariance of the KSS polynomials is 

coY{Xd{x),Xd{y)) = = [l + xy^. 

This fact suggest to homogenise the polynomials, that is, to introduce an auxiliary variable xq and to 
consider the polynomials 

d 

Xdixo,x) = 

n=0 

with xq, X G M. The polynomial X^ is homogeneous, that is X^(Xxo, Xx) = X^Xj^(xo, x) for any 
A; xo, X G M. Therefore, we can think of Xj^ as acting on the unit circumference S^. The covariance 
in this case gives 

cov(Xj^(xo,x),Xj^(yo,y}) = Y (xyr(xoyo)'^~'" = (xy + xcyo^ = ((xq, x), (yo,y}}‘^, 

where (,) stands for the usual inner product in 

Furthermore, denoting by Xy(I) the number of zeros of the process V on the set /; it is easy to 
see that 

2Nd(R) = Nj,o(S^). 

a 

Now, it is convenient to write X^ again as a polynomial in one variable by identifying (xq, x) G 
with the pair (sin(f), cos(f)) for some real t: 


d 


^d{t) := ^a„cos''(f)sin'^ "(f), 

n=0 


with real t. It is easy to see that 


Nd(R) = A^y-^([0, tt]) almost surely. 

In fact, X is a real root of Xd if and only if the radial projections of (1, x) onto S^, once identified with 
a point (sin(f), cos(f)), are roots of Vd; one of these projections correspond to f G [0,7r]. 
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Direct computations show that Yd is a centered stationary Gaussian process and that, for s, t G 
[0,7r], its covariance function is given by 

co\{Yd{s),Yd{t)) = ((cos(s), sin(s)), (cos(t), sin(t)))'^ = cos'^(t - s). 


Furthermore, 

coviYd{s), Yd{t)) = —dcos'^~^{t — s) sin(t — s), 

and 

COY{Yd{s),Yd{t)) = —d{d — 1) cos'^“^(t — s) sin^(t — s) + dcos^{t — s). 
In particular, Yar:{Yd{t)) = 1 and Var(yj(t)) = d for all t. 


Expectation of the number of roots of Yd'. 

Since the process Yd is stationary and for each fixed t the random variables Yd{t) and are 

independent centered Gaussian with variances 1 and d respectively, using Rice formula [3], we have 


E(iVy,[0,7r]) 


71 - ■PYa{0)W 



Vd. 


Time scale and covariance limit: 

The next step is to scale the time in order to get a limit behaviour for the covariances. It is convenient 
to use the unit speed parametrisation, so we define 


Zd{t) := Yd 



The number of real roots of Xd coincides almost surely with that of Zd in [0, x/dvr], that is 

Nd(J^) = Nz^{[0,VdTT]) almost surely. 

From now on, we restrict the process Zd to the interval [0, -v/dvr]. Since Yd is stationary, so is Zd- Let 
us denote by : [—\/d7r, ^/d-Tr] —>• M the covariance function of Zd, that is, rd{t) = cov {Zd{0), Zd{t)). 
It follows that 


rdit) = cos'^ 




( 1 ) 


Remark 1. Note thatvd is an even function and fort G [0,\/d7r/2] we have rd{y/d'K — t) = {—lYrd{t). 
This will imply that it suffiees to deal with Vd restrieted to the interval [0, \/d7r/2], as we shall see in 
Lemma 1. 


2 Asymptotic variance 


We need to prepare some preliminaries. 
Lemma 1. We have 


2\/d r'ZdTr/2 

nNd{Nd - 1)) = ^ / gd{t) 
^ do 


pl{t) + pd{t) arctan 


Pd{t) 


.\I^-Pdit)^ 


dt + 1. (2) 
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Here gd{t) = 27rpd{t)vd{t), being pd{t) = Pz^(o),Zd(t)(0, 0) the joint density of Zd{0), Zd{t) evaluated at 
(0,0); Vd{t) the conditional variance of Z'^{0) (and ofZ'^ft)) conditioned to Zd{H) = Zd{t) = 0 and pd{t) 
the conditional correlation between the derivatives Z'j^{0) and of Z'^{t) conditioned to Zd{Hj = Zd{t) = 0. 
We have 


Vdif) = 1 - 
Pd{t) = 


d cos^'^ ^ ( 4=1 sin^ ( 4^ 


\ \/d 


y/d 


1 (;fe) 

1 


27r Jl - cos2'i 


Pd{t) = cos 


d-2 


1- d sin^ ( 4^1 - cos^'^ ( 4^1 


\VdJ 


\VdJ 


V^/ 1 — cos^'^ ( 4^1 — dcos^‘^~‘^ ( 4^1 sin^ ( 4^1 


y/d 


y/d 


\ y/d 


Now, we pass to the asymptotic variance of Nd = Nz^{[0, -v/dvr]). We need the following asymptotics 
and bounds. 

Lemma 2. For each fixed t £ R, we have 


cos 




) d^oo 


Besides, these convergences are uniform in compacts. Furthermore, for 0 < a < 1 we have the 
following upper bounds 


cos'^ 



< 


g-atV2; 

cos'^(a); 


if 0 < t < ay/d, 
if ay/d <t< ■Ky/d/2. 


with a = 1 — a?/?> € (2/3,1). 

Remark 2. It is worth to say that this limit covariance defines a centered stationary Gaussian process 
X on [0,oo). The asymptotic behaviour of the number of real roots of Xd is intimately related to the 
asymptotic behaviour of the number of roots of X in increasing intervals. Similar situations occur in 
[1] (where this fact is indeed used explicitly to obtain the CLT) and [2j. The fact that the limit process 
X has Gaussian covariance function, and thus Gaussian spectral density, is remarkable. 

Nevertheless, we do not need this faet in the sequel. 

Lemma 3. For fixed t. 


gd{t) = 2 'k • Vd{t) • pd{t) 


1 - (1 + 

d—>oo (1 — g — i ^) 3/2 


=: 5(i); 


Pd/) ® 

d-^OQ 


-t2/2_ 


1 — — e 




1 — e 


=: p{t). 


Besides, 0 < g{t) < 1, \p{t)\ < 1 and g{t) -^t^o 0. Furthermore, there exists an integrable upper bound 
for the r.h.s. of Equation 2. 

Proposition 1 (limit variance for Nd). As d ^ oo we have 


Var(jVd) ^ ^2 2 F 

y/d d-^oo vr Jq 


ait) 


— p^(f) + pit) arctan 


pit) 


\/l 


-l\dt + l, 


where g and p are defined in Lemma 3. Furthermore, < oo. 
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Remark 3. Using Mehler formula we can write also 


a^ = 


poo 

I £ 


oo 2 

^2i M 




P {t){9{t) - Soi)dt, 


being 021 = 2(—1)^"’'^/(\/^2^£!(2^ — 1)) and 6 Kronecker’s delta. 

Proof. Recall that 

Var(iVrf) = E(iVrf(iVd - 1)) - {E{N)f + E{N). 
From Lemma 1 the normalized second factorial moment is 


EjNdjNd - 1)) 2 

Vd Jo 


\fd'K 12 


, - 2 rVd'^l'i ( p u\ 

gd{t)\/l - pi{t)dt + - gd{t)pd{t) arctan 


dt. (3) 


Let us look at the second integral in the r.h.s. of Equation (3). Let aVd < t < VdTTf2. We go 
back in our scaling: s e-)- t/y/d; so s € [a, 7r/2]. By the proof of Lemma 3 we have \pd{t)\ < cos'^“^(o) 
(see Equation (8) below for the details). Also by Lemma 3 gd{t) is bounded by constant. Hence, 


/ 


y/d-K 12 

y/n 


gdit)\pdit)\ arctan 


Pdit) 


'1 - Plit) 


\dt<[ 
I J a 


y/d-K l2 
■y/d 


cos*^ '^{a)dt -^d-^00 0 . 


Let 0 < t < ay/d. Lemma 3 gives the point-wise limit and the domination in order to obtain 


2 ^ ^ / Pdit) 

- / gd{t)pd{t) arctan 


2 f°° 

dt^d— / 5(t)/9(t) arctan 

^ Jo 


P{t) 


y/^-p"^/) 


dt 


In particular, this integral is finite. 

The first integral in the r.h.s. of Equation (3) cancel at infinity with {E{Nd)/. In fact, we can 
write 

2 i‘Udnj2 

(ENd/ = d = —y/d / dt. 

Jo 

Hence, the first integral in the r.h.s. of Equation (3) minus {E{N(i)/ gives 


2 

vr 


p'Ky/dj2 

! -£— 

0 r'Ky/dl2 

Vd{t)yJ 1 - pl{t) 

/ 

gd{t)yll- plit) - 1 


/ „ ^ 

Jo 


^ Jo 



dt. 


(4) 


From Lemma 3 it follows that the integrand of the l.h.s of Equation (4) tends to g{s) — 1. 

In order to obtain a domination, by standard manipulation, it follows that the important point is 
to bound the difference v‘^{t){l —p‘^{t))— l + r‘^{t) = — I + r‘^{t)) — {v‘^{t)p‘^{t)) in the numerator. 

The second term is easily bounded. For the first one, we have 


vj/) - 1 + rj{t) = 


1 - 


(icos2'^-2 ( 4. ) sin" ( 4. 


y/d 


y/d 


1 - cos^^ (^) 


-1-1- cos 


2d 


t 


After expanding the squares and cancelling the ones, we divide again into the cases 0 < t < ay/d 
and ay/d < t < y/dirfl. In the latter, the bound cos{t/y/d) < cos(a) suffices to obtain that the 
integral tends to 0. In the former, an uniform integrable upper bound for this difference follows 
easily by the triangle inequality and using the bounds d^ sin^^ {t/y/d) < t^^] cos{t/y/d) < e and 

cos’ ^{t/y/d) < cos ^(a). 
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The denominator is handled using the bound in Lemma 2. This gives the domination, so we can 
pass to the limit inside the integral. Hence, the first integral in the r.h.s. of Equation (3) minus 
(E(A^rf))^ tend, as d —>• oo, to 

2 

- / iait) - i)dt. 

^ Jo 

In particular, this integral is finite, thus, so is 

Finally, let us say that the convergence at 0 of the integral in Equation (4) follows from the proof 
of Lemma 3. The result follows. □ 

3 CLT 

Proposition 2. The normalised number of zeros of KSS polynomials converge in distribution towards 
a centered Gaussian random variable with variance cr^. 

The idea of the proof is the following: to write the normalised number of zeros as a chaotic series 
[8] and then to use the Fourth Moment Theorem [13] combined with Kratz-Leon method [9] in order 
to obtain the asymptotic normality. 

More precisely, we take the Ito-Wiener expansion of the normalised number of zeros [8] . Then, by 
Kratz-Leon method [9], the finiteness of the variance of allows to truncate the expansion and to 
derive its asymptotic normality from that of the sum of the first, say Q, terms. Finally, the Fourth 
Moment Theorem [13] gives a criterion to prove the asymptotic normality of the finite partial sums of 
the expansion. 

Proof. We apply Kratz-Leon expansion [8] to the processes on the interval [0,\/d7r]. Hence 

Nd - E{Nd) _ 

(^ 1/4 

q=2 

where 

i-y/d'K [ 9 / 2 ] 

h4 = -fT/A fq{Zd{t),Z'd{t))dt, fq{x,y) = '^bg_2ea2eHq_2e{x)H2eiy), ( 5 ) 

e=o 

with 021 = — 1)), bk = ^ip{0)Hk{0). Note that we can delete the term corre¬ 

sponding to O' = 1 since Lli(O) = 0; this is why we restrict our attention to zeros. 

We can express Ig^d as multiple stochastic integrals w.r.t B. In the first place, using a standard 
B.m. B we can write Zd{t) = f^hd{t, X)dB{X) with 

hd{t, /(O (^) (^) 

Then, from Equation (6), using the properties of the chaos and the stochastic Fubini theorem, see [2, 
Remark 2], we have Iq = Ig{gq{Xq)) with 

1 i"/d-K L 9 / 2 J 

= fjTj4 bq-2ja2j{h^'^~^^{s, Xq-2j) ® h'^^^s, X 2 j))ds; 

“ do j=0 

where and (8) stands for tensorial product. 

Now, to get the asymptotic normality of the standardised zeros, by Kratz-Leon method and the 
Fourth Moment Theorem, it suffices to prove that the contractions gq (8)fc gq{X 2 q- 2 k) tend to 0 in 
as d —)■ 00 for g > 2 and k = 1, - ■ ■ , g — 1 and X 2 q- 2 k G 
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Let Zk = (^ 1 ,..., Zk) and \ 2 q- 2 k = \-k ® ^q-k- The contractions are defined [13] as 

9q ^k 9qi^2q—2k) ~ I 9qi^kj ^q—k)9qi^k^ ^q—k)^^k- 
Actually, by the properties of stochastic integrals, we have Iq{gq{\q)) = {Sym{gq{\q))) being 

Sym{gq){\q) = 4 X] 9q{K), 

<T&Sq 

being Sq the group of permutations of {1,... , q} and A^- = ..., Xai^q))- So we compute contrac¬ 

tions for Sym{gq){Xq) instead of gq{Xq). 

Writing down the norm of the contractions is quite tedious, but the basic fact is that the isometric 
property of stochastic integrals implies that = r^{t — {t). Similarly, 

when the identified variable in the contraction involves the derivatives of h the resnlt involves the 
derivatives of r^. Taking this into account, it follows that 


\\Sym{gq) Sym{gq){X2q-2k)\\l 


1 

d 


ffff 

J J J J [O^y/d'K]^ 


E 

0<j<[q/2] 



E 

cr&Sq 




i=0 


s)T'{rd\t' 


s'))^^{r^d\s — s')y*{r'^\t — ds dt ds' dt'; 


where j = (ji, j 2 , js, j 4 ), vector inequalities are understood component-wise; Cj = 0^=1 
W = A = I3i{a,j), 7i = 7i(cr,j) and 6i = Si{a,j)-, ELi W = ELi A = k and Ei=i 7 * = 

Ei=i = q — k. Actually, there are some constrains for a,/3,j,6 with respect to j, (namely ai < 
(q — 2ji) /\{q — 2 ^ 2 ), Q !2 < (<? — 2ji) A 2j2 + (<? — 2j2) A 2ji, etc), but they are irrelevant for our purposes. 

We bound the covariances by their absolute value. Since Var(Zd(t)) = Var(Z^(t)) = 1, by Cauchy- 
Schwarz, each factor |r^*^(-)| < 1. Furthermore, since A: > 1 and <7 — A: > 1, we can bound from above 
the product of each group of factors (i.e.: with the same argument) by one of them. 

Hence, for some ii,i 2 ,h:H £ {0,1,2} we have 


\\Sym{gq) (8)fc Sym{gq){Xd)\\l < 


C 

7 


ffff 

J J J J [0,\/d7r]'^ 






“ '®0 ~ ^01 ds dt ds' dt', 


where (7 is a meaningless constant. Now, we make the change of variables: (x, y, u, t') 1 —)• (t — s,t' — 
s',s — s',t') and enlarge the domain of integration in order to have a rectangular one. Thus 


\\Sym{gq)^kSym{gq){Xd)\\l<— dt' \r^'^\x)\dx \rd^\y)\dy \r^^'^\u)\ du. 

Let us look at the three inner integrals. Note that since is even, so is the absolute value of 
its derivatives, so it suffices to integrate on [0, v 77 r]. Besides, since for t S [0, \/d 7 r/ 2 ] we have 
rd{\fdT^ — x) = {—l)'^rd{x), it follows that |r^*^(v 77 r — x)| = |r^*^(x)|,f G {0,1,2}. Then, we can 
further restrict the domain of integration to [0, v 77 r/ 2 ]. The finiteness of the integral then follows 
from Equation (1), by bounding the covariance by a polynomial (of degree at most 2) times cos'^{■/^/d) 
and then using Lemma 2. Hence, the contractions tend to 0. 

The result follows. □ 


Corollary 2. The asymptotic variance is strictly positive. 

Proof. From the Ito-Wiener expansion (5) it follows that = E ^2 T^ar(/q) with Var(/q) = lim^ Var(/q^rf). 
Thus, it suffices to prove that Var(/ 2 ) > 0. This is done exactly as in [3, Eq 10.42-10.43]. □ 
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Proof of Theorem 1. Put together Propositions 1 and 2 and Corollary 2. The approximated value for 
is obtained numerically from the formula in Proposition 1. □ 


It worth to say that this value is confirmed by simulations. 


4 Proofs of the lemmas 


Proof of Lemma 1. We compute the second factorial moment via Rice formula. By Equation (10.7.5) 
of [7] we have 


lE(iVd(A^d ^ I - i)9dit) [ \/l - P^it) + Pd{t) arctan | 




dt 


Denote 


f{t) = gd{t) [ ^Jl- pjit) + pd{t) arctan 


Pd{t) 




Then, f{y/d'K — t) = f{t). This follows from the properties cos(t) = cos(—t) = — cos(7r — t) and 
sin(t) = — sin(—t) = sin(7r — t) and from the fact that we only use even powers of the cosines and sines 
in gd and and that the signs also cancel in the product p arctan 

Then, using the change of variables x = y/dir — t in the interval [\/d7r/2, y/dir], we have 

s/dn 


ry a7T 

H^diNd - 1)) = 2 / /(t)(Vd7r - t)dt 

Jo 


= 2 
("s/dit 12 


p y/diT / 2 

/ f{t){\fd'i: — t)dt + 2 I f{t){y/d'K — t)dt 

Jo J y/div j 2 


('y/d'K 


rO 


= 2 / f{t){'/d7r — t)dt + 2 / f{Vd7r — x){'/d7r — ['/dTT — x]){—dx) 

Jo J y/div j 2 

r\/d7r/2 py/d-K l2 p\/d7rf2 


= 2 


py d'K j 2 pydTT/Z pydTT/Z 

/ f{t){'/dTT — t)dt + 2 / f{x)xdx = 2 / f{t)VdTTdt 
Jo Jo Jo 


□ 

Proof of Lemma 2. We start with the cosine. 

Assume that 0 < t < ay/d. Using Taylor-Lagrange expansion up to the first order for the cosine, 
we can write 


cos 


t 

y/d 


x\ 




= (1 + ^j ; withx = -y ( 1- 


sin(t*)t 

3y/d 


where t* S [0,t/y/d] C [0,a]. 

There exists c(a) such that \x/d\ < c{a) < 1. (in fact, if t < ay/d, then x/d = —a^/2 + sin(t*)a®/6. 
Since the summands have different signs, it follows that —1 < —a/j^ < xjd < a^/6 < 1. The claim follows with 
c(a) = max{a^/2, a'^/6} = a^/2. ) 

Again by Taylor-Lagrange expansion, for the logarithm this time, it follows that 

with X** G [0,x/d] C [0,a]. Hence, 










Note that log"(l + x**) = —(1 + x**) ^ G (—(1 — a) —(1 + a) ^]; so log"(l + x**) < 0. 

Now, for fixed t is easy to see that when d — >• oo then |t| < a\/d; X —t^l2 and log"(l + 
x**)x^I2d — 0. Thus, cos'^ (^) ^ as d —)• oo as claimed. 

It is also easy to check the uniformity of the convergence in compacts [0, (3] from the preceding 
expression. In fact, let t G [0,/3], then 



= 

gSin(t*)i3/6d3/2+log"(l+3;“(t))x2(t)/2d _ ^ 

KVd) 




< 


-t^l2 


max{ 


,sin(/3*)/33/6d3/2 _ 


,-(l-/32/2d)-2/34/4rf 


1 } ~^d 0; 


where we use the fact that the summands in the exponent have different signs. 


Now, we turn to the bounds. Since log"(l + x**) < 0 we see that cos^{t/y/d) < . Furthermore, 

since sin{-) < •, we have 0 < sin(t*)t/3\/d < a^/3. Hence, 


cos 


t 

y/d 


< e* = exp < — 


1 - 


sin(t*)t 

3y/d 


< exp <j — — ( 1 — 




(7) 


as claimed. 

Assume now that ay/d < t < yfdTrj2. In this case, we have cos{t/y/d) < cos(a) < 1, hence 
cos^{t/y/d) < cos'^(a). The result follows. □ 


Proof of Lemma 3. The limits are direct consequences of Lemma 2. 

Bounds; Let us look at the domination in a neighbourhood of t = 0. First, note that Pd{t) < 1 
since it is a correlation; so the sum is finite. 

Consider the factor 


2vr • Vdit) ■ Pdit) = 


1 — cos^'^ ( 4=1 — d cos^'^ ^ f 4^ 


\ y/d 


\ y/d 


Sin 


t 

y/d 


1_C0S2<^(^)) 


3/2 


By Taylor-Lagrange expansion, as in the proof of Lemma 2, we can write 


1 — cos 




. u 


+ tt + e“ — = -u - e“ —; 


where 


u = 2x + log"(l + ^** € [0)<7 [0,x]; 

^ = -^+ ^ G[0,t/A/d]. 

First, we bound x; since t lies on a neighbourhood of 0 (assume, for simplicity, that t < 1), the 
signature of the sinus is the same as that of t. Hence, 

-< X =-^- < -h —; <-h —. 

2 - 2 Qy/d ~ 2 6d - 2 6 

Now, we pass to the second term in the definition of u: since x < 0 and thus x** < 0, we have 

1 


-log"(l + x”) = 


(1 + X**/ 


> 1 . 
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Hence, 

_ lr,o-"n T-**'! , , , 

d d d \ 2 6 

On the other hand, x** > x > —1^/2, so 1 + x** > 1 — t^/2. Hence, 


- log"(l + X**)— > -r > 3 -^ + 77 • 


-log (l+g-)-<- 

Therefore, putting these bounds together 

+4 


- 1 '■2/o^2 - t 


{r/2y < 


-y- 




Hence, 


Finally, using that e“* < 1: 


1 


4(1 -tV2)2 


-y- 


<u<-r + 


1 / 

CS- 

to 

3 d [ 

2 


4(1 -t2/2)2' 


2 +4\ 2 

+ 


4(1 -t2/2)2 


<u< -r + —. 


i - 77 - 77 M + 


3 2 


4(1 -t2/2)2 


Thus 




2 2 

rt 9J 

< —tt-< 1 — cos^“ 


< 1 - cos^'^ ( ^ ) < -h 




4(1 -t2/2)2 

Similarly, we have cos^‘^~‘^ {t / \/d) > cos^^{t/^fd). Besides, 

t sin(r**)t2\2 


4(1 -t2/2)2- 


sin^(t/-\/d) = ( —= + 


\fd 


2d 


t*** G [0,t/Vd]. 


Hence, using the bound for 1 — cos2'^, we have 


- d cos^'^ ^ 


sin 


VdJ \Vd 


< —dcos 


2d 


< diy — 1 + 


t \( t sin(^***)^2^ 2 

y/d) \Vd 2d 

sin2(r**)t4 


4(l-t2/2)2y \d d3/2 


= - f + y 


t-l + 


4(1 -t2/2)2 


sin(r”) sin2(t*”)ti 


di/2 


4d 


4d2 


+ t + 


4(1 -t2/2)2 


< + 


4(1 -t2/2)2 


; if t < 1. 


Therefore, taking t < 1 
2vr • Vd{t) • pd{t) < 


f2 I 


+ {-t'^+t^+ 4 ( i 42 / 2)2 


t2 


V “I (^ + 


i3 

4(l-t^/2)^ 


3/2 


I J_ -L+t 

4(l-t^/2)^ 


4(l-t^/2)^ 


( 1 - ^ - 1 ( t 


3 2 1'' “I" 4(l-t^/2)^ 


3/2 


1 _ li _ 1 O j_ 

3 2 14(l-t^/2)^ 


3/2 ■ 


This gives an integrable (at 0) upper bound. 

Besides, Lemma 2 shows that the convergences are uniform in compacts. 
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Finally, it rests to obtain an integrable upper bound for large t (but recall that t G [0, \/d7r/2]). 
Assume that t > to- For t < a^/d we can use the bound cos'^(t/\/d) < Then, the upper bound 

is easy to obtain: 


27r • Vd{t) -pdit) < 


\Pd{t)\ < 


1 + (1 + 

(l_e-«F)3/2 ’ 

g-atV2 I + + e-"*" 

cos2(a) 1 - ‘ 


For t > ay/d is similar. Let s = 
we have 


t/y/d. Let us start with the correlation pd{t). Since cos(s) < cos(a), 


—dcos^'^ ^(s)sin^(s) + 1 — cos^'^(s) > —dco/'^ ^(a)sin^(a) + 1 — co/‘^{a) -^d 1; 
so, the demominator in pd{t) is positive for d large enough. Thus 


\pd{t)\ = co/ ^(s) 


|1 — (isim(s) — cos^“(s)| 

1 — cos^'^(s) — (icos^'^“^(s) sin^(s) 

1 — dsii/{s) — cos^'^(s) 


< cos*^ ^(s) 


= cos'^ ^(s) < cos'^ ^(o). (8) 


d-2/ 


1 — (isin^(s) — cos^'^(s) 


Similarly, since cos(s) < cos (a) in this region, we know that the conditional variance and the density 
(in this region) are bounded by constants (for instance, 1 and (27ry^l — cos^‘^(a))“^ respectively). 
This gives the necessary domination for the cases and concludes the proof of the Lemma. □ 
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